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Microscopic pyramidal pits in a reflective surface, a geometry similar to a retroreflector, are 
frequently used to enhance signal strength. The enhancement effect is generally attributed to surface 
plasmons, however, the sub-wavelength to near- wavelength dimensions of the pyramidal 3D geometry 
suggest contributions from diffraction and near-field effects. Our theoretical analysis of the light 
intensity distribution in the similar (but simpler) 2D geometry assuming a perfect conductor screen, 
that is, in the absence of any plasmon effects, shows that interference patterns forming within the 
cavity cause a significant resonant increase in local intensity. Such effect can be important for many 
applications, especially for the widely used Raman spectroscopy. Resonant enhancement without 
plasmons of the emitted Raman signal due to enhanced local field amplitude is also possible, which 
implies that the geometry practically implements a Raman laser. Comparison of diffraction patterns 
obtained with near-field and far-field approaches reveals that the near-field component is responsible 
for the observed dramatic intensity enhancement, and thus the Raman enhancement as well. 



PACS numbers: 42.25.Bs, 78.67.-n, 42.55.Ye, 42.65.Dr 



INTRODUCTION 

Micro- and nano-structured surfaces are frequently 
used to implement Surface Enhanced Raman Scattering 
(SERS), where interaction of the analyte with the sub- 
strate, mostly coupling between surface plasmons of the 
substrate and Raman modes of molecules 0, 0] , lead to 
the multiplication of the Raman signal intensity. Con- 
trolling the SERS enhancement, however, has proven 
problematic due to the random plasmon distribution. 

Nanovoids such as spheric cavities or grooves are 
known to host well-defined plasmon modes that inter- 
act strongly with incident laser light @, S| • It was thus 
proposed that nanovoids might be used to achieve tun- 
able Raman enhancement Q • One of the geometries that 
performs particularly well consists of an array of 1 /im 
deep pyramidal pits of ~90°opening angle jH]. Remark- 
ably, this mesoscale geometry has reflecting surfaces of 
a distance that scales between a few wavelengths and 
sub-wavelength and thus resembles an (unstable) optical 
resonator. This raises the possibility that in this par- 
ticular case the enhanced signal is generated by energy 
coupling into the Raman active medium from a high en- 
ergy standing wave, practically implementing a Raman 
laser. 

While parallel mirror micro-resonators have been stud- 
ied before 0, the tilted mirror geometry of the pyramidal 
pit requires a different formalism to account for diffrac- 
tion and near-field effects. Importantly, this geometry 
might be seen as a special case of an aperture (with an 
infinitesimally small opening at the tip of the pyramid). 
Light propagation through sub-wavelength apertures of 
various geometries is intensively studied due to the spe- 



cial properties of the optical near field, in particular, light 
localization and the accompanying intensity enhance- 
ment 171 la. Id. 1 1 Oil - While experimental evidence prooves 
its existence [lfj], the origin of the intensity enhancement 
is ambiguous as some theories attribute it to the excita- 
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tion of surface plasmon polaritons 
ers to intracavity (waveguide) resonances \\A 
other phenomena 17 , Qjl, Qjl] . While in the case of en- 
tirely sub-wavelength structures the comparable spatial 
distribution of plasmon- and diffraction-originated inten- 
sities prohibits an easy distinction, the dimensions of the 
pyramidal pit permit an attempt at the separation of the 
plasmon and diffraction effects. 

When assuming that the aperture is in a perfect con- 
ductor screen, in which the plasmon frequency is infinity, 
purely diffractive intensity distributions might be calcu- 
lated. The formulation we used is based on the near field 
treatment by Neerhoff and Mur [7j, l20j.We intend to find 
if resonant enhancement does exist in a pyramid-type ge- 
ometry, therefore we approximate the pyramidal pit with 
a V-slit, thus reducing the dimension of the problem to a 
plane, a simplification that, at the moment, is necessary 
for the solution of the problem. Importantly, the V-slit 
geometry is known to support standing waves as it can 
act as a waveguide; this application was discussed in a 
number of studies. [2l|, [H, Hi 0, fH 



THEORY 

In this section we summarize briefly the formulation 
used. For further details on a similar model see 0, S|- 
The configuration is shown in Fig. [TJ The slit is parallel 
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III. 



The following boundary conditions are imposed on G: 



dzGi(x,z)\ z _ b+ 
d z G 3 (x,z)\ z ^ - 



d x G 2 (x, z)\ 



x — >(^a,\ 



d x G 2 (x, z)\ 



for |a:| < oo, (2a) 

for |a;] < oo, (2b) 

for < z < b, (2c) 

for < z < b. (2d) 



Figure 1: Aperture model of the pyramidal pit. 



G\ and G3 must also satisfy the Sommerfeld radiation 
condition. Using these assumptions and conditions we 
may evaluate the three Green's functions. In regions I 
and III we find 



to the y axis in a screen of thickness b. The entrance 
width of the slit is a2, the exit width is a\ (in our case, 
di — > 0). The geometry is divided into three regions: I) 
|x| < 00, \y\ < 00, b < z < 00, II) \x\ < ai + a2 ^ ai z, 

\V\ < OO; 
Z < 0. 



< z < b and III) \x\ < 00, \y\ < 00, —00 < 



An incident plane wave propagating in the x — z plane 
in region I arrives at the slit at an angle 9 with respect to 
the z axis. The time harmonic magnetic field is constant 
and at the same time polarized in the y direction: 



H(a;, y, z, t) = U{x, z) exp(—iivt)e y . 



(1) 



In this case the electric field can be found from Maxwell's 
equations, and the diffraction problem is reduced to 
one involving the single scalar field U(x, z) in Eq. (Q]) 
in only two dimensions. In the three regions the field 
is represented by Uj(x,z) (j = 1,2,3). Due to time 
harmonicity the field satisfies the Helmholtz equation. 
In region I the field is composed of three components: 
Ui(x,z) = U l (x,z) + U r (x,z) + U d (x,z). Here U l (x,z) 
is the incident field which is assumed to be a plane wave 
of unit amplitude: U' l (x,z) = cxp[iki(xsm8 — zcosO)], 
U r (x,z) is the field that would be reflected without a 
screen: U r (x,z) = U l (x,2b — z). Finally, U d (x,z) is the 
diffracted field produced by the presence of the slit. 

For a perfectly conducting screen, the tangential E 
vanishes at the surface {d n U — 0). Further boundary 
conditions result from assuming that the diffracted fields 
and their derivatives vanish at infinity. The final bound- 
ary condition is the continuity of the field and its deriva- 
tive at the upper and lower slit boundaries. To find the 
field, the 2D Green's theorem is applied with one func- 
tion given by U(x, z) and the other by a standard Green's 
function. We assumed that the tangential E vanishes at 
the screen; therefore more is known about d n ll at the 
boundaries than is known about U itself. The Green 
functions for the V slit geometry can be found in |261 |. 



Gi{x, z; x 1 , z') = '-[H^ih^ + H^ihR')}, (3a) 
G 3 (x,z;x',z') = l -[H^(k 3 R)+H^(k 3 R'% (3b) 



where 



R = y/{x - X') 2 + (z - Z' f 



(4a) 

R' = ^{x-x') 2 + {z + z' -2b) 2 , (4b) 

R" = ^{x-x') 2 + {z + z') 2 . (4c) 

In region II, the method of images can be used: 



G 2 = 



OO 

i {H^[k\r-rl\] + H^[k\r-rl\]), (5) 



2 = 7TI 



where r' n and r" define a series of images reflected by the 
walls. The projection of an arbitrary point (x r , z') on the 
wall: 



z' + x' tana - 



1 



tana 



(x y 



tana 



(6a) 
(6b) 



From this the coordinates of the mirrored point are xm = 
2x x — x' and zm = 2z x — z' . A series of images of the 
original point might be generated with the rotation of 
the points by 4na (n being the number of rotation) . 

By using the definitions of the Green functions in the 
three regimes and the boundary conditions, one can find 
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the functions U d , U2, and U3: 



U d (x,z) 



U 3 (x,z) 



— Gi{x,z;x ,b)DUb(x )dx 

£2 
-a 2 

for b < z < co, 

ai 



63 
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Gz(x, z; x , 0)DUo{x')dx' 



U 2 (. 



x,z = - 



for — 00 < z < 0, 

G 2 (a;,«;»',0)£)f7o(a;') 
Uo(x')d z 'G2{x, z; x' , z')\ z ,^ +]dx' 

a 2 

- / [G? a (x,2;a5 / ,6)i3l7(,(a; / ) 



(7a) 



(7b) 



—Ub(x')d z rG2(x, z;x, z')\ z ,_ b -]dx' 
\x\ < ai + ai - ai z and — 00 < z < 0, (7c) 
where the boundary fields are the same as in 0]: 

U (x) = U 2 (x,z)\ z ^ +, (8a) 

DUq(x) = d z U 2 (x,z)\ z _ 0+ , (8b) 

U b {x) = U 2 (x,z)\ z ^ b -, (8c) 

DU b (x) = d z U 2 {x,z)\ z ^ h -. (8d) 

By applying the continuity equations to these functions, 
a set of four integral equations can be found. We solved 
the equations for U d , U2, and U3 numerically. 

RESULTS AND DISCUSSION 

We analysed the light diffraction in region II for sev- 
eral laser wavelength in the range that is frequently used 
for Raman spectroscopy. A typical result is shown on 
Figure [2] for A = 534 nm. Similar resonant intensity 
enhancement was observed for all wavelengths. This en- 
hancement is due to redistribution of electric charges (i.e. 
not plasmons) on the screen surface leading to standing 
wave formation with enhanced amplitude. The presence 
of the standing waves confirms that the geometry acts as 
an optical resonator. 

Another interesting feature is the formation of vortices 
in the cavity (Figure [3]). These vortices appear approxi- 
mately at z « nX/2, where n = 1, 2, 3, . . . It is important 
to note that at the surface the electric field vectors are 
mostly parallel to the walls. Thus, when considering the 
interaction with the Raman active material, the enhance- 
ment of different modes should occur in comparison to 
plasmon-enhancement where the electric field is mostly 
perpendicular to the surface. 

Far-field theory also predicts the standing waves and 
the vortices (Fig. [4]), however, a few remarkable differ- 
ences justify the use of near-field theory. First, at the 




-0.5 0.0 0.5 1.0 
x(@n) 

Figure 2: E x and E z electric and H y magnetic field dis- 
tributions, and the electric field vector in the V-groove for 
A = 534 nm. 
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Figure 3: Electric field vectors for three wavelengths: 434 nm, 
534 nm and 634 nm. 



slit entrance (z — > b~) a strong intensity suppression can 
be observed (see Figs. [5] and [3]) that cannot be seen in 
the case of a macroscopic slit with the same geometry 
(see Fig. [4]). Notably, far-field diffraction theory used to 
calculate intensity distributions in macroscopic objects 
does not contain the diffractive component that causes 
the suppression. Second, the maximum electric field oc- 
curs at the apex of the V-shape that confirms the obser- 
vations of Suzuki et al fl^ |. 

A comparison of the energy flux: x and z components 
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Figure 4: E x and E z electric and H y magnetic field distribu- 
tions (far-field theory) A = 534 nm. 



of the Poynting vector are depicted in Figure It is 
important to note that, while the standing waves appear 
at distances that satisfy the n\/2 criterium for the con- 
secutive wavelengths, the values of intensity maxima of 
the standing waves are strongly wavelength dependent. 
This dispersion effect might influence the efficiency of 
Raman-enhancement in these structures when using dif- 
ferent laser wavelengths, and could be used for the ex- 
perimental verification of the predictions of the calcula- 
tions. Note that enhancement at different wavelengths 
is not the same for electric field components and Poynt- 
ing vector components, e.g. Poynting vector shows no- 
table enhancement for A = 534 nm in comparison with 
A = 434 nm or A = 634 nm while for the same wave- 
length the maximal value of the electric field vector is 
the smallest. 

In summary, theoretical analysis of light diffraction in a 
sub-wavelength to near-wavelength V-groove in a perfect 
conductor screen was performed by using near-field ap- 
proach. It was shown that the geometry acts as an optical 
resonator where the intensity enhancement is determined 
by the near-field effects, i.e. without plasmons. The field 
vectors form vortices around the intensity maxima, being 
mostly parallel to the walls of the groove, which deter- 
mines the Raman modes of surface bound molecules en- 
hanced in this configuration. The intensity enhancement 
is wavelength dependent, thus different laser wavelengths 
might lead to different Raman enhancement when using 
this geometry. Comparison of results obtained with near- 
field and far-field approaches reveals that the near-field 



component is responsible for the observed dramatic in- 
tensity enhancement, and thus the Raman enhancement 
as well. 
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Figure 5: Poynting vector for three wavelengths: 434 nm, 534 nm and 634 nm. Note the "moving" of the standing wave. 
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